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A numerical and systematic parameter study of three-dimensional vesicle electrohydrodynamics 
is presented to investigate the effects of different fluid and membrane properties. The dynamics 
of vesicles in the presence of DC electric fields is considered, both in the presence and absence of 
linear shear flow. For suspended vesicles it is shown that the conductivity ratio and viscosity ratio 
between the interior and exterior fluids, as well as the vesicle membrane capacitance, substantially 
affect the minimum electric field strength required to induce a full Prolate-Oblate-Prolate transition. 
In addition, there exists a critical electric field strength above which a vesicle will no longer tumble 
when exposed to linear shear flow. 


I. INTRODUCTION 

The interaction of lipid vesicles with external electric 
fields has been a subject of growing interest from both 
an application and a mechanistic standpoint W- One 
of the most striking phenomena in vesicle electrohydro¬ 
dynamics is the Prolate-Oblate-Prolate (POP) transition 
mm- Consider a prolate vesicle where the major axis is 
initially aligned with the applied electric field. If the fluid 
encapsulated by the vesicle has a lower conductivity than 
the surrounding fluid there is a chance that the vesicle 
will transition to an oblate shape, where the major axes 
are perpendicular to the electric field. Over time the vesi¬ 
cle will evolve back into the prolate shape (aligned with 
the electric field) and stays in that equilibrium condition 
permanently 00- In order to observe the full POP tran¬ 
sition the external electric field has to be strong enough 
to overcome the forces of the vesicle membrane and the 
hydrodynamic forces of the fluid. 

In physical experiments capturing a full POP transi¬ 
tion with a single DC pulse is not a trivial task P3 |H]. 
Strong fields can easily cause membrane poration which 
causes the vesicle to stay in the oblate shape permanently 
or to collapse before reaching the final prolate condition. 
Therefore, applying an electric field with the appropriate 
strength and duration is crucial in directing the behavior 
of the vesicle into the desired dynamics and configura¬ 
tion. 

Although a considerable number of experimental works 
in the past decade have been reported on vesicles exposed 
to electric fields a thorough numerical analysis 

of vesicle electrohydrodynamics is quite scarce. Existing 
analytical studies are solely limited to nearly spherical 
vesicles 0 0 or planar membranes m nu and make 
the Stokes assumption for the fluid flow. In this paper 
a systematic parameter study is carried out for a three- 
dimensional vesicle in the presence of DC electric fields. 
A full Navier- Stokes solver is employed which is able to 
capture possible inertia effects in the vesicle dynamics 
arising from the higher Reynolds number of electrohy¬ 
drodynamic problems compared to typical hydrodynamic 
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simulations in the absence of electric fields. The effects 
of fluid and vesicles properties on the electrohydrody¬ 
namic behavior are examined and the critical parameter 
thresholds required to observe a full POP transition are 
determined. This includes the effect of membrane capac¬ 
itance and of the viscosity and conductivity ratios of the 
embedded and surrounding fluids for a range of different 
reduced volumes. Vesicles in the presence of combined 
shear flow and DC electric fields are also investigated and 
the critical field strength required to change a normally 
tumbling vesicle into a tank-treading one is determined 
for different viscosity ratios. 

II. THE MATHEMATICAL MODEL AND 
NUMERICAL METHODS 

The model used in this paper is based on the com¬ 
putational approach presented by the authors in Ref. 
HU- Consider a three-dimensional vesicle suspended in 
an aqueous solution where the inner and outer fluids may 
have differing fluid viscosity and conductivity. The mem¬ 
brane is assumed to be inextensible with constant en¬ 
closed volume and surface area. The sum of the mem¬ 
brane forces is given as 

-Kcdfi-2 HK + V 2 a H)n + 7 Hn-V 3 'Y, (1) 

where n is the outward facing normal vector on the inter¬ 
face, 7 is the interface tension, k c is the bending rigidity, 
H is the total curvature (sum of principle curvatures) 
and K is the Gaussian curvature at the interface. For a 
vesicle in the presence of an external electric field these 
forces are balanced by the sum of the hydrodynamic and 
electric stresses. The hydrodynamic stress at the inter¬ 
face is given as 

T hd = n • [Thd] = 71 • {T^d — T h d) > (2) 

where T^d = — pi + /i(Vu + V T u ) is the bulk hydrody¬ 
namic stress tensor while /i, u and p refer to viscosity, 
fluid velocity, and pressure, respectively. The + sign is 
used when the interface is approached from the outer 
fluid while the — sign represents the interface being ap¬ 
proached from the inner fluid. Square brackets, [ ], indi¬ 
cate the jump of a quantity across the interface. 
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The electric field stress at the interface is expressed as 

T el = n ' [T el] = Tl ‘ (T+i — T el ), ( 3 ) 

with the Maxwell tensor, T e j, given as T e i = e(EE — 
\ (E • E) J), where E = — V^> is the electric field, e the 
permittivity, and is the electric potential field. 

By employing a level set method to implicitly track 
the interface varying fluid properties at any point in the 
domain can be written in a single relation, e.g. p(x) = 
/i _ + (/i + —/i _ )H(0(cc)) where T-L is the Heaviside function 
and (p is the level set function. One is therefore able 
to write the momentum equations of binary fluids with 
varying properties into one single formulation. 

Introduce to, uo, and as the characteristic time, 
velocity, length, and electric field. Three important di¬ 
mensionless quantities can then be defined as Ca = 
/i + (l + r])a 3 / (t 0 ft c ), Mn = t 0 e + ^o/(^ + (l + v)) and 
Re = pu^to/fi^. The quantity Ca indicates the strength 
of the bending forces, while Mn indicates the strength 
of the electric field effects and Re is the Reynolds num¬ 
ber. This allows for the single-fluid formulation of the 
Navier-Stokes equation to be written as [T2] 


Du 

*T>t 


- Vp + • (/i (Via + V T u)) 

+ 8 (</>) II V0|| (V s 7 - 

+ chu 5 w(?r~ 2KH+ * H )*4' 

e (^EE -^(E-E)I^j 


+ ^^)||V^|| 


n. 

(4) 


Additionally, the fluid incompressibility and surface in¬ 
extensibility constraints are given as V • u = 0 and 
V s • u = 0, respectively. Note that in Eq. all quan¬ 
tities are nondimensionalized and the singular contribu¬ 
tions of the bending, tension, and electric field forces have 
been transformed into localized body force terms using 
the level set and Dirac delta functions. All fluid prop¬ 
erties are normalized with respect to the external fluid 
values. 

The leaky-dielectric model is used to obtain the electric 
potential in the domain. The electric potential is the so¬ 
lution to the Laplace equation in each region, V 2 ^ >± = 0, 
with an electric potential jump at the interface due to 
the capacitive property of the membrane. This poten¬ 
tial jump is called trans-membrane potential (V m ) and 
evolves over time by the non-dimensional relation, 

8V 

C m -^+u-W (C m V m ) = n-E+, (5) 

where C m = (C m a)/(t o<s + ) is the non-dimensional mem¬ 
brane capacitance, C m is the actual capacitance of the 
membrane and A = s~/s + is the conductivity ratio be¬ 
tween the two fluids. Note that any trans-membrane con¬ 
ductance is currently ignored. In this work an immersed 
interface method is used to take into account the time- 
varying trans-membrane potential and the discontinuous 


fluid conductivities. This is done by deriving jump con¬ 
ditions for the electric potential and its first and second 
derivatives in order to achieve second order accuracy. For 
details of the numerical implementation interested read¬ 
ers are referred to Ref. m ■ 

For the electrohydrodynamic computations with no 
imposed shear flow the characteristic time, to, is set 
to the membrane charging time scale i.e. to = t m = 
(a + W In simulations with combined shear flow 
and an external electric field the time scale associated 
with the applied shear flow is used i.e. to = tj = 
where 70 is the applied shear rate. Typical experimen¬ 
tal values reported for physical properties are given as 
a ~ 20 p m, k c ~ 10 -19 J, s + ~ 10 -3 S/m,, e + = e~ « 
10 -9 F/m, Cm « 10 -2 F/m 2 , p « 10 3 kg/m 3 , y 0 « 1 
s _1 , p~ = /i + « 10 -3 Pa s and for the strong electric 
field Eo = 10 5 V/m 0 HU H5] . Using these values the 
dimensionless parameters in the presence of an external 
electric field and in absence of shearing flow are found to 
be Re = 0.19, Ca = 3.8 x 10 4 , Mn = 18 and C m = 0.1. If 
the vesicle under the combined effects of electric field and 
imposed shear flow is considered, then the corresponding 
dimensionless parameters are Re = 10 -3 , Ca = 10 and 
C m = 2 x 10 -4 . Note that unlike the cases with 70 = 0 
where strong electric fields are applied, in the combina¬ 
tion of the two effects only weak electric fields are consid¬ 
ered therefore, the dimensionless strength of the electric 
field is within the range of Mn = 10. 

The moving interface is modeled through the use of 
a semi-implicit, gradient-augmented level set jet scheme 
which is an extension of the original jet scheme intro¬ 
duced in D3I- The semi-implicit algorithm is imple¬ 
mented in three steps. The first step is the Lagrangian 
advection of the cell centers. In the second step the val¬ 
ues at the cell centers are smoothed using a semi-implicit 
scheme. This is done based on introducing a smoothing 
operator similar to the idea used in Refs. mm- In the 
third step the center values are projected back to grid 
points and the jet of the solution including the gradients 
are updated; see Ref. [12] for implementation details. 

To tackle the hydrodynamic problem a recently devel¬ 
oped Navier-Stokes projection method is used. Both lo¬ 
cal and global fluid incompressibility and surface inexten¬ 
sibility conditions are satisfied in the numerical method 
by splitting the pressure and tension correction terms 
into spatially-constant and spatially-varying components 
and solving for the system of unknowns through the use 
of a Schur complement decomposition strategy. For de¬ 
tails of the numerical approach and the special discretiza¬ 
tion see Ref. na- 


in. RESULTS 

In all the present simulations the vesicle surface area 
is fixed to 47 r and the initial shape is a prolate ellipsoidal 
vesicle. The computational domain is a box with the size 
of [—4.5,4.5] 3 . Periodic boundary conditions are taken in 
the x— and z—directions while wall-boundary-conditions 
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are taken in the y— direction. The electric field is im¬ 
posed by setting an electric potential equal to <f> = — y 
on the wall-boundaries. In the case where shear flow is 
also applied the velocity at the wall-boundaries is taken 
to be u = y. Note that the electric field and fluid velocity 
have been normalized by Eq and uq. A Cartesian collo¬ 
cated mesh with uniform grid spacing is used and the 
grid spacing is chosen to be h = 0.075 with a minimum 
time-step of At = 5 x 10 -4 . Further information about 
the domain size, time step, and grid spacing choices can 
be found in Ref. [12]. 

The vesicle configuration is described by the deforma¬ 
tion parameter: 


where a is the axis length of the vesicle in the direction 
parallel to the applied electric field and b is the length of 
the axes perpendicular to the electric field, see Fig. [1] A 
prolate shape is given by D > 1 (major axis aligned with 
the direction of the electric field) while the oblate shape 
has D < 1 (major axis perpendicular to the direction 
of the electric field). The deformation, D , is directly 
calculated from the interface location as described by the 
level set. 


strengths (represented by Mn values) are presented in 
Fig. H In the simulation results in the bottom series a 
strong electric field of Mn = 28 is applied while a weaker 
field, Mn = 20, is used for the results in the top row. 
In the former case, the electric forces are strong enough 
to overcome the membrane forces of the vesicle and the 
viscous forces of the fluid and this drives the vesicle into 
the oblate state where D < 1. However, when the weak 
electric field is applied the initial prolate vesicle transi¬ 
tions into to a slightly less prolate profile but is unable to 
obtain the prolate to oblate transition. The final shape 
is the same in both cases. This indicates that there ex¬ 
ists a critical electric field strength, Mn c , above which a 
full POP transition occurs to the suspended vesicle ex¬ 
posed to an external DC field. As will be shown below, 
this critical electric field strength depends on a number 
of material parameters. 
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FIG. 1, Sketch of the prolate and oblate shapes and their 
orientations with respect to the external electric field. 


FIG. 2. Simulation results showing the full POP transition 
(bottom series) vs. a vesicle that remains in the prolate 
regime (top series). The reduced volume in both simula¬ 
tions is v = 0.93 and other simulation parameters are set 
to Re — 0.19, Ca — 3.8 x 10 4 , Cm — 0.1, A — 0.1 and 
rj = 1.0. The colors on the surface (online version) indicate 
the distribution of the trans-membrane potential over time, 
with darkest blue (lower pole) indicating a membrane poten¬ 
tial of Vm = — 1.6 and darkest red (upper pole) indicating a 
membrane potential of V m = 1.6 for the saturated membrane 
capacitance. 


The electrical properties of the two fluids, including 
the permittivity and electric conductivity, are among the 
most important factors in determining the dynamics of 
vesicles exposed to external DC electric fields. The ratio 
of the electric permittivity to the electric conductivity in 
each fluid gives the bulk charge relaxation time expressed 
as t± = This time represents the rate at which each 
fluid can supply charges to the interface which is different 
for fluids with different electrical properties. The ratio 
of the charge relaxation time of the interior fluid to the 
exterior fluid is given by t~/t+ = £/A where ( = e - /e + 
and A = s~/s + . This measure is of major importance in 
vesicle electrohydrodynamics and shows the competition 
between the two fluids in conducting charges towards the 
membrane. For typical vesicle experiments it is reason¬ 
able to assume £ = 1 and hence the deviation of the 
conductivity ratio from one is the factor that determines 
different vesicle dynamics. 

To illustrate the types of deformations in the full POP 
transition and how it compares to an incomplete transi¬ 
tion the shape of a vesicle with two different electric field 


To study the role of the fluid conductivity ratio in POP 
dynamics the critical Mason number, Mn c , is determined 
for various conductivity ratios. This behavior is illus¬ 
trated in Fig. [3] for a range of different reduced volumes. 
The region above each curve shows where the POP tran¬ 
sition occurs while the region below the curves indicates 
the Mn values that are not strong enough to force the 
vesicle to full POP transition. As the the conductivity 
ratio A increases, larger and larger electric field strengths 
are required to induce the POP transition. Theoretically 
speaking, for a vesicle with A = 1, i.e. matched con¬ 
ductivity between the two fluids, the Mn c approaches 
infinity. In fact, as will be demonstrated in Sec. [m]it 
is not possible for a vesicle with match conductivity to 
achieve a POP transition. 

The membrane capacitance is another important fac¬ 
tor in determining the dynamics of the vesicle in a DC 
field. In particular, the capacitance appears as a primary 
parameter in the time scale associated with the mem¬ 
brane charging given, t m . The larger the capacitance the 
larger duration of time electric field forces can act on 
the membrane. In the context of the POP transition an 
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FIG. 3. The effect of the conductivity ratio, A, on the POP 
transition. The region above the curve shows where the POP 
transition happens while the region below the curve repre¬ 
sents a vesicle that remains prolate for all the time. The 
transition is plotted for four different reduced volumes, v = 
0.90,0.93,0.95, and 0.98. The two regions are separated by 
the critical field strength, Mn c where the dynamics changes 
from purely prolate to prolate-oblate-prolate. The nondimen- 
sional parameters are set to Re = 0.19, Ca = 3.8 x 10 4 , 
rj — 1.0 and Cm = 0 . 1 . 

interesting comparison can be made between the mem¬ 
brane charging timescale and the electrohydrodynamic 
time scale of the applied field given by t e hd = fI € +~e^ • 
The electrohydrodynamic time scale largely depends on 
the electric field strength: a strong electric field is associ¬ 
ated with a shorter electrohydrodynamic time. In theory, 
for a nearly spherical vesicle a full POP transition will 
occur if t e hd < t m [Z] - In such situations electric field 
forces have sufficient time to induce a full POP deforma¬ 
tion before the trans-membrane potential is saturated at 
time t m . As the capacitance increases the charging time 
of the membrane capacitance, £ m , increases as well which 
implies that even with a weaker electric field (larger t e hd) 
a full POP transition is still possible as long as t e hd is suf¬ 
ficiently smaller than tm. 

Note that the membrane capacitance affects the 
amount of deformation as well. An increase of t m due 
to larger capacitance leads to larger deformations as the 
electric forces only influence the dynamics of the mem¬ 
brane in proportion to the time that they are allowed to 
act. Since the capacitance is a major factor in determin¬ 
ing this time the implication is that the deformations are 
smaller for smaller dimensionless capacitance, C m , and 
deformations grow larger as C m does. This manifests 
the importance of the membrane capacitance in the type 
and magnitude of deformations. 

Figure [4] shows how changes in the membrane capac¬ 
itance leads to a vastly different critical field strengths, 
Mn c , needed to observe a full POP transition. This be¬ 
havior is consistent across all investigated vesicle reduced 


volumes. For the four reduced volumes investigated the 
critical field strength scales as the inverse of the mem¬ 
brane capacitance, Mn c oc C” 1 . This relationship will 
be further discussed in Sec. EO 



FIG. 4. The effect of the membrane capacitance, C m , on 
the POP transition. The region above the curve shows where 
the POP transition occurs while the region below the curve 
represents a vesicle that remains prolate for all the time. The 
transition is plotted for four different reduced volumes, v = 
0.90, 0.93, 0.95, and 0.98. The two regions are separated by 
the critical held strength, Mn c , where the dynamics change 
from purely prolate to prolate-oblate-prolate. The nondimen- 
sional parameters are set to Re = 0.19, Ca = 3.8 x 10 4 , 
7] — 1.0 and A — 0.1. 

The effect of the viscosity ratio (77 = /i - //i + ) on the dy¬ 
namics of the vesicle is shown in Fig. [5] for four different 
reduced volumes. Numerical experiments demonstrate 
that as 77 increases a stronger electric held is needed to 
overcome the resisting hydrodynamic forces of the fluids. 
A larger viscosity ratio is associated with a larger viscous 
damping force which in turn leads to slower dynamics and 
thus a larger Mn c required to induce the full POP tran¬ 
sition. The critical electric held strength appears to be 
linearly related to the viscosity, Mn c oc 77 . 

One consequence of these results is that it appears that 
given the same electrical properties, smaller reduced vol¬ 
umes require a stronger held to achieve a full POP tran¬ 
sition. This is most likely dependent on the initial shape, 
as smaller reduced volumes begin as more prolate (more 
stretched out) vesicles. This means that as the reduced 
volume decreases larger deformations are required to ob¬ 
tain the prolate-oblate transition. Future work will inves¬ 
tigate this, as it is possible to have many different initial 
vesicle shape conditions for the same reduced volume. 

Next, the combined effects of an electric held and shear 
how are considered. In this case the strength of the ex¬ 
ternal electric held is increased from zero to investigate 
the changes in the dynamics of the vesicle. Note that 
here the shear how time scale is used as the characteristic 
time-scale. With typical experimental values this results 
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FIG. 5. The viscosity ratio effect on the vesicle behavior tran¬ 
sition from the prolate regime (region below the curve) to 
POP regime (region above the curve). The two regions are 
separated by the critical field strength, Mn c where the dy¬ 
namics changes from one to another. This transition is plotted 
for four different reduced volumes, v = 0.90, 0.93, 0.95 and 
0.98. The non-dimensional parameters are set to Re = 0.19, 
Ca = 3.8 x 10 4 , Cm = 0.1 and A = 0.1. 

in non-dimensional parameters Re = 0.001, Ca = 10 and 
Cm = 2 x 10 -4 . For typical situations the shear flow 
time scale is approximately = Is while the membrane 
charging time scale is approximately 10 -3 s [12]. There¬ 
fore, in this situation the trans-membrane potential is 
taken to respond instantly to changes in the membrane 
configuration and is thus at a pseudo-steady-state. 

The effect of the electric field is to align the vesicle 
vertically with the direction of the field. Therefore, for 
a vesicle in the tumbling regime an increase in Mn will 
result in a slower tumbling period. A close observation 
of vesicle dynamics under the combined shear and an ex¬ 
ternal weak electric field reveals that the periodic behav¬ 
ior in this situation is more akin to vacillating-breathing 
motion where the vesicle does not follow a rigid-body¬ 
like rotation anymore, but instead the poles retract and 
additional deformations are induced in the shear plane. 
If the field is strong enough, the electric field forces act¬ 
ing against the shearing forces of the fluid can eventually 
bring the vesicle into the tank-treading regime. 

These different dynamic regimes are shown in Fig. 
[ 6 ] In this example the vesicle has a reduced volume of 
v = 0.93 and the viscosity ratio is r] = 10. This ratio 
is larger than the critical viscosity ratio needed to ob¬ 
tain tumbling in shear flow when in the absence of an 
electric field. The top row in the figure shows the dy¬ 
namics over time when the electric field is not present 
in the simulation (Mn = 0). As would be expected the 
vesicle tumbles with the profile of the vesicle remaining 
nearly constant and rotating periodically in time. The 
vacillating-breathing-like motion is seen for the results 
shown in the middle row. In this simulation a weak elec- 



FIG. 6 . Dynamics of a vesicle under the combined effects 
of shearing flow and an applied electric field at various times. 
The reduced volume is v = 0.93 and a viscosity ratio of 77 = 10 
is used. Other simulation parameters are Re = 10 -3 , Ca m 10 
and Cm = 2 x 10 -4 . From left to right the snapshots are 
given at t — 0.0, 2.5,6.0, 7.0,8.0,10.0,15.0,17.0. The top row 
shows the tumbling vesicle in the absence of external electric 
held (Mn = 0). The middle row corresponds to a vesicle 
with a vacillating-breathing-like motion under an electric held 
strength of Mn = 3. The bottom row shows a tank-treading 
vesicle under an strong electric held strength of Mn =15. 

trie held of Mn = 3 is used. Unlike the tumbling motion, 
in this regime the prohle of the vesicle undergoes changes 
over time as it periodically rotates about the vorticity 
axis. Finally the bottom row demonstrates the vesicle 
dynamics of the same vesicle under an electric held with 
a stronger strength (Mn = 15). The vesicle reaches an 
equilibrium, tank-treading inclination angle and stays in 
that position permanently. 

The critical Mason number at which the transi¬ 
tion from tumbling-like motion (either tumbling or 
vacillating-breathing) to tank-treading occurs is denoted 
as Mn 7 . The required Mn 7 for two different reduced 
volumes over a range of viscosity ratios is plotted in Fig. 

m 


IV. DISCUSSION 

To better understand why a vesicle will undergo a 
prolate-oblate-prolate transition in the presence of an 
electric held let us consider the electric held forces acting 
on the membrane. For simplicity and to match the com¬ 
putational results assume that the electric permittivity 
ratio is one, e + = e ~. Next, consider only the normal 
electric held forces on the interface. Tangential electric 
held forces will drag the fluid into motion but do not 
induced changes in the vesicle membrane. 

It is straightforward to write the normal electric force 
in terms of the normal and tangential components of the 
electric held, 

r e m=C ((. E + 2 - E- 2 ) - (. E+ 2 + E + 2 -Ef - Ef )), 

(7) 

where E t and E 5 are the tangential and E n is the normal 
electric held at that point. Note that the notation here 
is E+ 2 = (. E+f and E~ 2 = (. E~f. 

This clearly shows that both normal and tangential 
electric held components contribution to membrane de- 
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FIG. 7. Tumbling-like to tank-treading transition for vesicles 
with two different reduced volumes, v = 0.90 and v = 0.93, 
under the combined effect of imposed shear flow and weak 
DC electric fields. The non-dimensional parameters are set 
to Re = 0 . 001 , Ca = 10 , x =*1, C m = 2 x 10“ 4 and A = 0 . 1 . 
Membranes with Mn ^ 0 have a pseudo-steady-state trans¬ 
membrane potential. The region above the curve is the tank- 
treading regime while the region below the curve shows the 
tumbling-like regime. Increasing the field strength from zero 
gives rise to a damping force against the motion of a normally 
tumbling vesicle. 


formation. To further simplify the discussion, consider 
the two points A and B located respectively at the top 
pole and the equator of the vesicle, see Fig. [ 8 j Due 
to their particular locations with regard to the external 
field the normal electric field force at these two points is 
simplified to 

r^n=y(£ n +2 -£- 2 ), (8) 

tT n = (J Ef + Ef - Ef - Ef). (9) 

One of the conditions on the interface is the continuity 
of the normal electric field, s + E+ — s~E~ = 0. It is 
therefore possible to write E+ = A E~ and thus at point 
A the normal electric force becomes 


-A 

T el 


^ + (A-1) r _2 


n = 


EZ 


( 10 ) 


After the electric field is applied charges from both 
the enclosed and surrounding fluids start migrating to¬ 
wards the interface. However, due to the differing electric 
conductivities the charges accumulate at differing rates. 
This leads to an apparent charge, 


Q = e+(E+-E~), ( 11 ) 


assuming matching permittivities. 




FIG. 8 . A sketch of the electric field, induced charges and 
electric forces when an initially prolate vesicle is transitioning 
to the oblate shape (left figure) and when the oblate vesicle 
is evolving back to the equilibrium prolate transition (right 
figure). 


Let the electric potential be activated at time 0. The 
electric fields inside and outside the vesicle are initially 
matched and non-zero: E+ = E~ ^ 0, E^~ = E^ 0, 
and E^ = E^ ^ 0. According to Eq. Q the force along 
the equator is zero. The direction of the force at the vesi¬ 
cle poles depends on the conductivity ratio. When A < 1 
the normal force is in the — n direction, which is a verti¬ 
cally compressive force. On the other hand, when A > 1 
the force is extensional and thus the POP transition can 
not occur. This demonstrates the necessary condition of 
A < 1 when the fluid permittivities are matched. 

Now consider what occurs when t oo. Given enough 
time charges are able to accumulate on both sides of the 
interface. Therefore Q = 0 and the electric field inside 
the vesicle becomes zero, E~ = 0 [7;. Accordingly, for 
matched permittivities E~ = E+ = 0 while E = 0 and 
E^ = 0. Note that in this case E+ / 0 and E£ / 0. 
From Eq. (§ this results in a zero normal force at point 
A while along the equator, point T>, the normal force 
reduces to a purely compressive force. Notice that this 
behavior is independent of the fluid electrical conduc¬ 
tivities, indicating that the final shape will always be a 
prolate one. 

Finally, a rational behind the influence of the three 
main parameters investigated here, membrane capaci¬ 
tance, fluid conductivity ratio, and fluid viscosity ratio, 
on the critical electric field strength needed to obtain a 
POP transition is presented. First examine the contribu¬ 
tion of the membrane capacitance. The membrane capac¬ 
itance influences how fast the trans-membrane potential 
reaches a steady state value, Eq. (§. The time it takes 
for the trans-membrane potential to reach steady state is 
inversely proportional to C m . As the amount of time that 
the electric forces can act on the membrane is directly re¬ 
lated to the membrane-charging time, a decrease in t m 
results in a system which does not have enough time to 
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respond and obtain a full POP transition. Therefore the 
Mason number must increase at a rate proportional to 
C~ l 2 3 4 to ensure that the POP transition occurs. 


Next consider the viscosity ratio, 77 . It was stated in 
Sec. [ml that for a POP transition to occur the elec¬ 
trohydrodynamic time, t e hd, needs to be smaller than 
the membrane charging time, t m . As the electrohydro¬ 
dynamic time increases linearly with an increase in the 
viscosity ratio, t e hd oc I + 77 , the membrane charging time 
must also increase linearly to r] to maintain the t e hd < t m 
relation. Recall that for results in the absence of shear 
flow the characteristic time was set to the membrane 
charging time, to = t m , and therefore the Mason number 
depends linearly on t m . Thus, an increase in the viscos¬ 
ity ratio results in a linear increase in t e hd, which results 
in a linear increase of t m so that the POP transition can 
occur, which finally leads to a linear increase of Mn c . 


The conductivity ratio presents a more complicated 
picture. Based on the previous discussion only A < 1 will 
result in a POP transition and thus all cases with A > 1 
are ignored. Examining t m it is clear that an increase 
in the conductivity ratio will result in a decrease of the 
membrane-charging time. This is due to the fact that 
as A —)> 1 the fluids are able to supply charges to the 
interface at the same rate and thus a charge mismatch 
does not occur. Following the discussion regarding the 
viscosity ratio, a corresponding decrease in t e hd would be 
required for a POP transition to occur. This is achieved 
by increasing the electric field strength and therefore in¬ 
creasing Mn. Analysis of the equations results in a scal¬ 
ing of Mn c oc (2A)/(2 +A), which can not account for the 
rapid increase of Mn c as A approached 1, Fig. [3J It is 
thus necessary to turn to the forces acting on the mem¬ 
brane and how they relate to A, Eq. (10). Clearly, as 
the conductivity ratio approaches 1 the forces acting on 
the membrane asymptotically approach zero. It is there¬ 
fore expected that rapid growth of the electric field would 
be needed obtain a full POP transition. It is suspected 
that at small conductivity ratios the increase in Mn c is 
dominated by (2A)/(2 + A) while as A —>> 1 the growth 
becomes dominated by (1 — A) -1 . Further investigations 
are needed to verify this. 


V. CONCLUDING REMARKS 

In this paper a parameter space study of 3D vesicles 
exposed to DC electric fields was presented. In partic¬ 
ular, the focus of this work was on the effects of differ¬ 
ent fluid and membrane parameters on the transition be¬ 
tween prolate and oblate shapes for a suspended vesicle, 
in addition to the transition from the tumbling to tank- 
treading regime of a vesicle under the combined effects 
of a DC electric field and shear flow. 

The electric field strength, indicated by the critical Ma¬ 
son number, Mn c , required for a POP transition was de¬ 
termined. A decrease in the reduced volume of the vesicle 
was always accompanied with an increase in the critical 
Mason number. While an increase in the viscosity ratio 
between the inner and outer fluids brought about an al¬ 
most linear increase of Mn c , the change in the dynamics 
of the vesicle by varying the conductivity ratio between 
the two fluids and the membrane capacitance of the vesi¬ 
cle was found to be more dramatic and substantial. 

The results of the vesicle behavior under the combined 
effects of shear flow and weak DC electric fields revealed 
the remarkable influence of the electric field in changing 
the standard behaviors of tank-treading and tumbling 
vesicles. Investigations showed that the application of 
the electric field generates a damping force against the 
motion of the normally tumbling vesicle. If this force 
is strong enough the tumbling vesicle stops the flipping 
tumbling behavior and undergoes a tank-treading motion 
The critical Mn 7 for this transition to occur was deter¬ 
mined over a range of different viscosity ratios for vesicles 
with two different reduced volumes. 

The results and analysis presented here assumed 
matched fluid permittivities and only considered a single 
initial vesicle shape. It is expected that different initial 
vesicle shapes could result in differing dynamics while 
the equilibrium shape would remain the same as pre¬ 
sented here. Additionally, a relaxation of the matched 
permittivity condition could result in behavior not ob¬ 
served here. Future work will explore the influence of 
initial vesicle shape and permittivity differences on the 
dynamics of vesicles exposed to DC electric fields. 
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